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HIAD & Flexible TPS Background 
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• Potential missions have been identified that will require a 
planetary entry system to have an aeroshell much larger in 
diameter than the diameter of any feasible launch vehicle. 

• The inflatable aeroshell of a HIAD entry system relies on a 
flexible thermal protection system (TPS) to prevent the 
underlying structure from exceeding its thermal limits. 

• Candidate flexible TPS materials and layups have undergone 
extensive aerothermal arc jet testing in a stagnation 
configuration to evaluate thermal performance and provide 
boundary condition and in-depth temperature measurement 
data for thermal model correlation and validation. 

• Outer fabric protects the underlying insulation layers from 
being directly exposed to the incident convective heat flux 
and the aerodynamic shear forces. 

• The insulator layers reduce thermal soak back. 

• Gas barrier prevents hot gas impingement on to the 
underlying structure. 
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• Given uncertainties in material properties, environments, etc. 
Prevent bondline over temperature of the F-TPS 


• Use probabilistic tools to enhance F-TPS 
performance and increase survivability 
during off nominal entry conditions, 
while reducing F-TPS margined mass 


Provide more knowledge and insight 
about the F-TPS design and its 
performance 



A! A A SciTech Forum January 08, 2015 


4 



• Uncertainties captured without stacking worst cases together 

• Decoupled approach two separate Monte Carlo simulations are run 

■ Trajectory 

> Output is a set of dispersed trajectories 

■ Thermal model 


> Each resulting trajectory from the trajectory Monte Carlo is used as input to the 
thermal model Monte Carlo 
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Monte Carlo Simulation - Implementation 
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• Margin Process 

■ Thermal Monte Carlo simulation run using several layups, each time 
increasing the number of insulation layers 

■ Peak bondline temperature tracked and presented as a bondline 
temperature distribution 

■ Margined thickness found by selecting the layup which satisfies mission 
reliability requirement for the F-TPS 



Layup #1 


Layup #2 
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Probabilistic Sizing Process Overview 





3<j = 375°C 3cr=21 6°C 
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Determine Bondline 
Temperature PDF 
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Bondline Temperature Distributions 
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Probabilistic Sizing Overview 


Sizing Objective 

• Bondline Temperature Limit: 400°C 



• Margin Policy Dictates: 97% chance Bondline temperature will be maintained below 
400°C (Predicted +2o bondline temperature must be less than 400°C) 


Bondline Temperature Prediction Requirements 


Aerothermal Environment 
(CFD) 


Trajectory Dispersions 
(POST) 


Material Thermal Response 
(COMSOL) 


Thickness Margin 

• Material Response Uncertainty - 

Parametric, Stochastic, & Structural 

• Trajectory Dispersions - heat flux and 
heat load 

• Aerothermal Uncertainty - heat flux 



Outer Fabric 



1 st Gen Flexible TPS 


Insulator 


Bondline 
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Sensitivity Analysis 
(Monte Carlo #1) 



Sensitivity Analysis 


Thermal Response Model Overview 




Stagnation Test 



Inside Model Holder 


Instrumented Flexible TPS Test Specimen 


Thermal Response Model Layout 
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Sensitivity Analysis 



Parameters Dispersed in Monte Carlo Sensitivity Analysis 


Parameter 

Distribution 

Type 

Mean or 
Nominal 
Value 

Standard 
Deviation 
or Range 

Insulator Virgin Conductivity (k vrg ) SF 1 
(Low Temperature Range) 

Normal 

1.0 

0.2 

Insulator Virgin Conductivity (k vrg ) SF 2 
(High Temperature Range) 

Normal 

1.0 

0.2 

Insulator Char Conductivity (k chr ) SF 1 
(LowTemperature Range) 

Normal 

1.0 

0.2 

Insulator Char Conductivity (k chr ) SF 2 
(High Temperature Range) 

Normal 

1.0 

0.2 

Insulator Specific Heat (C pc ) SF 

Normal 

1.0 

0.15 

Ap (Virgin Density - Char Density) SF 

Uniform 

1.0 

0.8 to 1.2 

Gas Specific Heat (C pa ) SF 

Uniform 

1.0 

0.1 to 1.0 

Gas Viscosity (p a ) SF 

Uniform 

1.0 

1.0 to 2.0 

Latent Heat of Reaction 1 (h R , ), SF 

Uniform 

1.0 

-30 to 30 

Insulator Porosity (cp ) SF 

Uniform 

1.0 

0.575 to 1.09 

Insulator Permeability (k Y ) SF 

Normal 

1.0 

0.2 

Outer Fabric Conductivity SF 

Normal 

1.0 

0.2 

Backing Material Conductivity SF 

Normal 

1.0 

0.2 
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Sensitivity Analysis 


Sensitivity Analysis Results 

Time Evolution of Input Parameter Contributions to Bondline Temperature Uncertainty 
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Uncertainty Breakdown: Temperature BL 
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Thermal Model Correlation 
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Temperature (degC) 



Thermal Model Correlation 


Correlated Model Predictions 


Correlation Parameters 


Arc Jet Test Results 


Imposed Temperature Boundary Conditions 


Constant Heat Flux 



Markers => Arc Jet Data 
Solid Lines => Model Predictions 



Parameter 

Scaled 

Value 

Insulator Virgin Conductivity (k wo ) SF 1 

0.885 

Insulator Virgin Conductivity (k wo ) SF 2 

1.09 

Insulator Char Conductivity (k phr ) SF 1 

0.97 

Insulator Char Conductivity (k rhr ) SF 2 

1.01 

Backing Material Conductivity SF 

1.75 

Outer Fabric Conductivity SF 

1.70 


Modeled Gen 1 Layup 

Measured Surface Pressure Boundary Condition 
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Bondline Temperature 


Time (sec) 
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Correlated Model Predictions 


Markers => Arc Jet Data 
Solid Lines => Model Predictions 


Constant Heat Flux 
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Thermal Model Correlation 


Correlated Model Predictions 


Markers => Arc Jet Data 
Solid Lines => Model Predictions 



Lower 

Constant Heat Flux 

2 Insulator Layers* 


Modeled Layup 


Measured Surface Pressure BC 

i 



TC-2R (BC) 
TC-3R 

TC-4K 

TC-5K 


TC2 

TC3 

TC4 


2 


Time (sec) 


Time (sec) 



Profile Heat Flux 


Modeled Layup 


Measured Surface Pressure BC 

l 



-TC-2R (BC) 
TC-3R 
TC-4K 
-TC-5K 

-TC-6K 

-TC-7K 





Uncorrelated 


Correlated 


A! A A SciTech Forum January 08, 2015 


16 



Error Analysis Results - 12 Arc Jet Cases Compared to Predictions (4 profile and 8 constant HF) 


Error Indicator 

Uncorrelated 

Correlated 

Average Max BL 
Temperature Difference °C 

165 

30 

Average Normalize BL 
Temperature RMS% 

27 

5 


Max Bondline Temperature Difference (Predicted - Measured) °C 

I H 4 •— H 4 1 1 Hi • 1 1 — +—-\ 

-50 -40 -30 -20 -10 0 10 20 30 40 50 
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Uncertainty Analysis and Sizing 
(Monte Carlo #2, #3, & #4) 
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Uncertainty Analysis & Sizing 



Monte Carlo Input Parameters and Uncertainty Quantification 


Parameter 

Distribution 

Type 

Mean or 
Nominal Value 

Standard 
Deviation (lo) 
or Range 

Insulator Virgin Conductivity (k wa ) SF 1 

Normal 

0.885 

0.106 

Insulator Virgin Conductivity (k wo ) SF 2 

Normal 

1.09 

0.106 

Insulator Char Conductivity (k rhr ) SF 1 

Normal 

0.97 

0.131 

Insulator Char Conductivity (k rhl .) SF 2 

Normal 

1.01 

0.131 

Insulator Specific Heat (Cps) SF 

Normal 

1.00 

0.0738 

Latent Heat of Reaction 1( h R1 ) SF 

Uniform 

1.00 

0 to 2 

Heat Flux (q nu; ) SF 

Normal 

1.00 

0.07 

Bondline Temperature Bias 

Uniform 

0°C 

0 to 50 

Insulator Thickness SF 

Normal 

1.00 

0.084 

Insulator Density SF 

Normal 

1.00 

0.083 
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Uncertainty Analysis & Sizing 



Monte Carlo Input Parameters: Trajectory Dispersions 

120 

100 
80 
60 
40 
20 
0 

0.96 0.97 0.98 0.99 1 1.01 1.02 1.03 1.04 1.05 

Normalized Heat Load Histogram 

• Two thousand (2000) off-nominal trajectories were provided from a separate Monte Carlo simulation run using POST 
where initial orbital conditions, vehicle aerodynamic parameters, atmospheric conditions, and other parameters 
affecting the reentry trajectory were dispersed. 

• The boundary condition parameters derived from each of the 2000 trajectory dispersions as inputs into each of the 
corresponding 2000 thermal model Monte Carlo samples. 
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Uncertainty Analysis & Sizing 


Sizing: Nominal Case, Four Insulator Layers 




Modeled Layup 


Heat Flux and Surface Pressure BC 

outer fabric 
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Maximum Bondline Temperature Histogram with Lognormal PDF Fit 
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Uncertainty Analysis & Sizing 


Sizing: Four and a Half Insulator Layers 




Modeled Layup 


Heat Flux and Surface Pressure BC 


outer fabric 
outer fabric 
insulation 


4 



insulation 


insulation 
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Maximum Bondline Temperature Histogram with Lognormal PDF Fit 
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Uncertainty Analysis & Sizing 


Sizing: Five Insulator Layers 



Modeled Layup 

Heat Flux and Surface Pressure BC 

outer fabric 

outer fabric 

insulation I I 

insulation I I 

insulation i i 

insulation I I 

insulation I I 

gas barrier 
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Uncertainty Analysis & Sizing 



Maximum Bondline Temperature Uncertainty Contribution Breakdown (5 Layer Case) 
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• The probabilistic sizing methodology prescribes one additional insulation layer on top 
of the nominal thickness to mitigate parametric, stochastic, and structural uncertainties. 

• A thickness margin of one layer of insulation was shown to result in a 99.87% probability 
of maintaining the bondline below 400°C, exceeding the given margin policy. 

• This thickness margin may shrink as knowledge is gained about the parameters 
which are contributing the significantly to the uncertainty in the maximum 
bondline temperature or by developing a higher fidelity thermal model, enabling a 
higher quality correlation and a lower maximum bondline temperature bias. 

• The probabilistic sizing process works well and it has enabled the opportunity to make 
informed TPS design and mission risk tradeoffs. 

• A HIAD project may decide to decrease the TPS thickness margin to one half a layer 
of insulation resulting in a 5% probability of bondline over-temperature as 
determined by the second Monte Carlo simulation for uncertainty analysis. 
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• Lesson Learned: Sensitivity analysis showed that the gas specific heat (C pg ) has a large 
influence on the maximum bondline temperature, this indicates that advection is a 
significant mode of heat transfer. 

• The specific heat of the gas, as well as other gas properties, can be better defined 
by introducing multi-species gas mass conservation into the thermal model and 
determining the properties of the mixture. 

• This will reduce the thermal response model structural uncertainty and allow for 
a better correlation to arc jet test data which will decrease the bondline 
temperature bias applied in the uncertainty analysis Monte Carlo simulations. 

• This probabilistic sizing process can be readily implemented on other TPS layup 
designs provided a validated thermal response model, appropriate boundary condition 
dispersions, and thermal model input parameter uncertainty quantification. 
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• Demonstrate this process using a higher fidelity model with multi- 
species-gas-mass conservation integrated in to the numerical 
solution approach. 


* Further investigate uncertainties and the correlations between the 
thermal response model input parameters which were dispersed 
or not dispersed in this study. 

* Employ automated correlation by inverse analysis (Grant 
Rossman work). 


• Implement this process for sizing the TPS of the Terrestrial HIAD 
Orbital Reentry (THOR) aeroshell (HIAD technology demonstration 
flight test). 
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Back Up Charts 
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Types of Uncertainties (1 of 3) 


• Parametric uncertainties 



■ Uncertainty that results from errors in the 
underlying physical models, or model input 
parameter estimates 

>i.e. char thermal conductivity, pyrolysis gas enthalpy, 
decomposition kinetics constants 

❖Properties that are difficult to measure are often estimated or 
theoretically calculated 

>This type of uncertainty can be reduced through 
testing and analysis but not eliminated 
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Types of Uncertainties (2 of 3) 
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• Stochastic variability 

■ Uncertainty that arises from natural fluctuations 
in the physical environment, or from material 
property variability 

>i.e. Varying atmospheric conditions, lot-to-lot property 
variation 


>This type of uncertainty is always present and can not 
be reduced 

❖Typically accounted for in the entry trajectory simulations 
using a Monte Carlo simulation 

■ These uncertainties are what most people think of when 
they think about uncertainties 
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Types of Uncertainties (3 of 3) 
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• Structural uncertainties 

■ Arise in numerical simulations due to the fact that 
simulations employ mathematical models to simulate 
physical phenomenon 

■ Mathematical models often make simplifying 
assumptions which are sometimes only valid over a 
limited range of conditions 

■ Computational truncation errors 

■ i.e. Tauber-Sutton radiative heating correlation, 
transfer coefficient-recovery enthalpy approach for 
convective heating boundary condition (FIAT/CMA), 
differential equation discretization method 
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* Biases are factors applied to the nominal thickness 
to account for known deficiencies in the 
methodology used to calculate the nominal TPS 
thickness 

■ Biases account for structural uncertainty and are applied to 
both nominal and margined thicknesses 

■ A bias should be applied if the thermal model shows a consistent 
difference between prediction and test data 

■ i.e. CEV TPS applied a bias to thickness based on the calculated 
recession 

■ For HIAD we need to add a bias to account for an uncorrelated 
model 

* Margin is applied to the nominal thickness to 
account for the parametric and stochastic 
uncertainty types 
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• Random variables are statistically correlated when the probability 
distribution of one depends on the value of another random variable 

• The strength of the dependence between two random variable is 
assessed by calculating their correlation coefficient which varies 
between -1 and 1. A values of 0 means no correlation. 

• In semi-plain English, when a functional relationship exists between 
two parameters, they are correlated and can not be treated as 
independent random variables 


Sample correlation coefficient 



Where x and y are the mean values 
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Random Variable Correlation 
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• This may seem like a problem, but its not 

• For each sample in Monte Carlo simulation, a vector of random 
numbers are generated corresponding to the thermal model 
parameters of interest, {R} 

• If there are parameter correlations present, {R} must be modified 
based on each parameters correlation with one another. The 
correlation coefficients are collected in a correlation matrix [C] 

• The correlation matrix [C] can’t be used directly, so an 
intermediate matrix [U] must be found such that: 

[Uf[U] = [C] 

• Once [U] is found by using Cholesky decomposition, the 
correlated vector of random variables can be determined by: 


{R C } = [U]{R} 
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Random Variable Correlation 
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• In reality there is a good chance that correlations exist between 
density and many other thermal model parameters, however . . . 

■ In the COMSOL model, we are not able to determine whether or not 

correlations exist between density and thermal conductivity, density and specific 
heat, density and permeability, etc, because there is no functional relationship 
established between these parameters in the COMSOL model 

• The correlations can be determined experimentally, but at an expense 

• By adding the degree of char to the COMSOL model, we may be able 
to computationally determine if there are correlations, but it is not as 
simple as just tracking the density change over temperature since 
gas conduction plays a significant role in the effective thermal 
conductivity of the material 
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• Currently TPS is sized using a root-sum-square 
(RSS) approach. Typically using three cases 

>Case 1 : size heat shield to nominal bondline temperature, 
apply trajectory dispersion multiplying factors to heat flux 

>Case 2: size heat shield to nominal bondline temperature, 
apply trajectory dispersion and aerothermal uncertainty 
multiplying factors to heat flux 

>Case 3: size heat shield to toiverthan nominal bondline 
temperature, apply trajectory dispersion multiplying factors 
to heat flux 

■ These three cases are RSS’ed together to arrive at 
the margined thickness 

T RSS = 1-10 * [Tj + V (t 2 - T i) 2 + (t 3 - T J 2 
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